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SUMMARY 


Measurements are presented of the aerodynamic forces and moments 
acting on a wing-tank configuration, with or without fins, oscillating 
in pitch about the wing-root midchord. The reduced-frequency range was 
from 0.050 to O. 657 , whereas the Mach number and Reynolds number ranges 

were from 0.18 to 0.75 &nd 0.9 X 10^ to 9»5 X 10^, respectively. 

Comparisons of the experimental aerodynamic forces and moments and 
their respective phase angles acting on the wing-tank combination with 
those acting on the wing alone indicated that the overall lifts and moments 
were greater when the tank was attached. The forces on the tank alone 
were compared with those determined by an engine -nacelle theory developed 
by Andropolus, Chee, and Targoff in Air Force Technical Report Number 6353. 
The agreement between measured and calculated lift was very good. 


INTRODUCTION 


As a result of the uncertainties and difficulties involved in the 
analytical treatment of the unsteady aerodynamic forces acting on bodies 
of revolution such as tip tanks, the validity of flutter analyses for 
configurations equipped with external stores has been subject to question. 
Flutter has occurred for wings equipped with tip tanks, and analysis has 
indicated the possible importance of the air forces contributed by the 
tanks in influencing the flutter characteristics. 


In order to further the knowledge of the unsteady aerodynamic forces 
on wing— tank configurations, measurements have been made of the air forces 
acting on a low-aspect-ratio wing equipped with a tip tank oscillating in 
pitch about the midchord. In order to determine the effect of tank fins, 
the measurements were made with two different fins and were compared with 
measurements made with the tank in a clean condition. For comparison, 
averaged experimental data from a previous investigation for a wing with 
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an aspect ratio of 2 are shown with the current results of the measure- 
ments of the wing— tip-tank combination. The coefficients for the oscil- 
lating tank in the presence of an oscillating wing are compared with cal- 
culated coefficients for an engine nacelle as determined from a paper by 
Andropoulos , Chee, and Targoff (ref. l). 

The aerodynamic coefficients are presented for a range of reduced 
frequency from 0.050 to 0.657 and for a range of Mach numbers from 0.18 

to 0 75 The Reynolds number range is from 0.9 X 10^ to 9*5 X 10 . These 
measurements were made in the Langley _ 2- by 4-foot flutter research tunnel 
by using a. resonant-oscillation technique* 


SYMBOLS 


A 

a 

a l 

B,C,R,Y,A 

c 

c r 

El 

g 

St 

§vac 

Ie 


i = 7- 1 

K 

K s 


aspect ratio 
lift-curve slope 

response amplitude to be determined 

dummy variables 

chord of wing, ft 

reference chord 

bending stiffness 

structural damping coefficient 

damping coefficient of wing in air stream 

damping coefficient of wing in a near vacuum 

effective inertia of oscillating system accounting for 

2 

dynamic deformation of system, ft-lb-sec 
correction 

effective spring constant of oscillating system, ft -lb /radian 


k 


reduced-frequency parameter, ux:/2V 
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L oscillating-lift vector acting on the wing for oscillations 

about the midchord axis, positive when acting upward, 

L = | l| ^(rt+Ax/180) . ^3^1 + ll2 ) e io* 

L(x) distribution of air load 

^1 nondimens ional coefficient of 

displacement 

2 2 nondimensional coefficient of 

velocity 

| 2| = Jl i 2 + 

M Mach number 

M a osc illat ing-moment vector acting on wing for oscillation 

about midchord, referred to axis of rotation, positive 
for leading edge up, 

M a = | M a | e 1 ( <wt+ ^2 It / 1 8o) = ^5 | | a | ^ + im^j e 1 ^ 

m mass of wing 

effective reference mass 
% mass at wing tip 

m l nondimensional moment coefficient in phase with angular 

displacement 

m 2 nondimensional moment coefficient in phase with angular 

velocity 

P applied aerodynamic loading 

p intensity of applied loading 

q dynamic pressure, lb/sq ft 

S area of wing, sq ft (taken as 1 square foot when used in 

coefficients ) 

t time, sec 

V velocity of test medium, fps 


lift in phase with angular 
lift in phase with angular 
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P 

T 

01 

02 

<h 

O) 

% 

^Vac 


coordinate of wing from root to tip 
amplitude of wing deflection 
unit tip amplitude 

angle of attack at midsemispan station as a function of time, 
|a|e ia) ^, radians 
root shear 

mass density of test medium, slugs/cu ft 
tip of wing 

phase angle that lift vector leads angle of attack, 
tan-1 (l 2 / *1), deg 

phase angle that pitching-moment vector leads angle of attack, 
tan'-^n^m^j , deg 

phase angle that load leads root shear 

circular frequency of oscillation of wing-tank combination, 
radians/ sec 

circular frequency of first natural wing bending, radians/sec 

circular frequency of oscillation of wing in a near vacuum, 
radians/ sec 


Subscripts : 

V? wing in presence of tank 
T tank in presence of wing 
WT wing-tank combination 


Dots indicate differentiation with respect. to time. 
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APPARATUS AND METHOD 


Tunnel 

The Langley 2- by 4 -foot flutter research tunnel was used for the 
present investigation. The test mediums were air or Freon -12 as noted 
in table I. The use of Freon -12 as a test medium permits the attainment 
of approximately twice the reduced frequency obtained in air for a given 
Mach number and frequency. 


Model 

The model consisted of a 12 -inch-span, 12 -inch-chord wing with a 
5 . 2 - inch-diameter tank, 24 inches long, mounted over the wing tip. (See 
figs. 1 and 2.) Fabricated construction was employed for the wing; a 
steel box spar carried four evenly spaced ribs to which plywood skin was 
attached in order to form an NACA 65AOIO airfoil section. The wing was 
designed to have high natural frequencies in order to reduce the amount 
of correction to the measured forces due to aeroelastic deformations and 
to bending inertia loads. (See appendix A.) The first natural cantilever 
bending frequency of the wing alone was about 130 cycles per second. 

The tank was constructed of thin balsa nose and tail sections attached 
to a thin magnesium center section with ordinates as shown in figure 2. 

The tank was mounted on a strain-gage balance attached to the wing tip and 
had a gap with the wing surface of about l/8 inch. This gap was necessary 
to isolate the tank from the wing so that only the tank forces and moments 
would be sensed by the strain-gage balance in the tank. The exposed area 
of the wing was 0.80 square foot and the plan-form area of the tank was 
0.69 square foot. 

In order to study the effects of tail fins attached to the tank, some 
tests were made with the tank equipped with one of two fins. One fin was 
trapezoidal in shape and the other fin was delta in shape. (See fig. 2 .) 
Both fins had the same area ( 0.08 square foot) and were mounted so that 
their centroids of area were the same distance rearward of the pitch axis. 

The frequency of the tank on its strain-gage beams was 117 cycles 
per second, the first bending frequency of the wing-tank combination was 
80 cycles per second, and the first torsion of the wing-tank combination 
was over 150 cycles per second. 

The semispan model was mounted as a cantilever beam in an oscillator 
mechanism at the tunnel wall. This mechanism permitted the model to oscil- 
late in pitch about the midchord axis. The model was mass-balanced about 
the axis of oscillation in such a way that there were no lift reactions 
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when the model was oscillated in a near vacuum. During the tests, the 
frequency of the pitching oscillation of the system varied from 19 to 
lJ+ cycles per second. 


Oscillator Mechanism 

The oscillator mechanism may be considered as a simple, torsional, 
vibratory system and is illustrated in figure 3. The system consists of 
a torsion spring which is fixed at one end with the other end attached to 
a bearing- supported, hollow steel tube to which the model and base plate 
are clamped. The mechanism was oscillated in torsion at the natural re- 
quency of the simple spring-inertia system by applying a harmonically 
varying torque through the shaker coils attached to the steel tube. In a 
near vacuum, this natural frequency was approximately 20 cycles per second. 

The bearings were contained in housings which were carried on columns. 
These columns were designed to include stress -sensitive regions and were 
equipped with strain gages from which the aerodynamic lift would be deter- 
mined The vertical reactions at the fixed end of the torsion spring were 
negligible because of the flexibility of the torque rod and the small 
deformations experienced by the strain-gage columns. 

The electromagnetic shakers consisted of stationary coils furnishing 
a steady magnetic field and moving coils which were attached to the steel 
tube The moving coils were driven by a variable-frequency oscillator. 
These moving coils were alined so that the direction of the applied forces 
was perpendicular to the direction of the lift. Tims, even if a pure 
torque were not applied to the steel tubes, no resulting force would be 
sensed by the lift gages. Provision was made for interrupting the power 
to the moving coils in order to obtain a power-off decaying oscillation. 


Ins t rument at ion 

The instrumentation was designed to provide continuous signals that 
were proportional to the lift and angular position and to provide a means 
of measuring their amplitude and time relationship. The lift reactions 
were converted to electrical signals by means of wire strain gages attached 
to the supporting columns. The gages were connected so that only lifting 
loads were sensed. An electrical signal from a wire strain gage mounted 
on the torsion spring so as to sense torsional strains was calita-ated to 
give the angular displacement in terms of the wing incidence. The t 
was attached at the wing tip through a strain-gage dynamometer so that 
the tank forces could be separated from the total forces on the combina- 
tion. The signals were filtered to eliminate noise and higher harmonics, 
and their magnitudes were measured by using a vacuum-tube voltmeter. 
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The angular -position signal was recorded on a recording oscillograph 
during the decay of the oscillation for the purpose of obtaining the 
damping coefficient. 

The phase measurements were made with an electronic counterchrono- 
graph. The time lapse between a given point on a lift (or tank moment) 
signal and a corresponding point on the position signal was measured while 
the wing-tank combination was oscillated at a constant frequency. The 
period of oscillation was measured with this instrument by determining 
the time lapse between corresponding points on the same signal. 


Calibration 

The angular position of the wing was dynamically calibrated with the 
signal from the torsion strain gages by a photographic technique. Time 
exposures were taken of a fine chordwise line on the outer edge of the 
tank for various amplitudes while the strain-gage output was recorded. 

The amplitude of oscillation of the wing was obtained from the envelope 
position of the line on the outer edge of the tank and correlated with 
the strain-gage signal. By using this procedure and a line on the leading 
edge of the wing, it was determined that, at the maximum frequency of 
oscillation (20 cycles per second), the tip angle of incidence exceeded 
the root angle of incidence by less than 1 percent. 

The signals from the balance columns were calibrated in terms of 
pounds of force per unit of signal strength. Known loads were applied 
to the wing, and the column reactions were determined by treating the 
wing shaft system as a simple beam with overhang. (See, for instance, 
fig. 3.) The reaction forces were then related to the respective signals. 
The vacuum-tube voltmeter was calibrated dynamically by using a low- 
frequency oscillator and a voltage divider. These dynamic calibrations 
were then related to forces by use of the open-circuit load calibration 
of the strain gages. The readings from the vacuum-tube voltmeter are 
believed to be within ±4 percent of true signal. 

The phase-measuring system was calibrated at various frequencies by 
using standard resistance-capacitance phase-shift circuits and by using 
a cam-operated set of cantilever beams on which strain gages had been 
mounted. The latter system had distortion and noise and approximated the 
worst tunnel condition. Calibrations of the phase meter indicated that 
the phase angle may be determined within ±3° of true value with a noisy 
signal and within t2° of true value with a clean signal. In order to 
minimize errors in phase introduced in the electrical operations, a tare 
value of phase was obtained at each reading by applying either the lift 
or angular-position signals through both channels of the electrical 
circuits . 
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Data Reduction 


The total lift forces as received from the balances contain an aero- 
dynamic component and an inertia component which arise from the bending 
deformation of the wing. In order to correct the measured lift to the 
aerodynamic lift, it was necessary to correct for the inertia forces due 
to wing deformation. A discussion of this correction is given in appen- 
dix A. The inclusion of this correction leads to a factor which, when 
multiplied by the measured lift, gives the actual applied lift. The value 
of this factor is O.98. In order to estimate possible error incurred by 
neglecting the aerodynamic forces and moments arising from the bending 
deformation, these forces were included in the analysis in appendix A and 
were found to be less than 1 percent of the correction due to wing defor- 
mation caused by the inertia forces. The phase angle contained a 

component due to these forces that tended to increase by less than 

1°. The moments due to the bending deflection were found to be negligible 

j_v 0 to the magnitude of the measured moment. Since the aerodynamic 
effects due to wing bending were within the accuracy of the measurement, 
no effort was made to adjust for these quantities. 

The in-phase moment on the wing-tank combination was determined from 
the change in resonant frequency due to air flowing over the wing and tank 
as indicated in reference 2. Since the torsional damping was small, its 
effect on the frequency is neglected and the entire shift is attributed to 
the in-phase moment. The moment coefficient in phase with the angular dis- 
placement is given by 


The dependency of m x on the small difference of two quantities of 
approximately the same magnitude leads to considerable loss in accuracy 
and consequent scatter in the data. 

The quadrature -moment coefficient for the wing-tank combination was 
determined by operating on the time history of the angular position 
obtained during the power-off decay of the oscillation. The moment coef- 
ficient in phase with the angular velocity is given by 




The derivation of this equation is treated in appendix B. 
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The phase angle 02, WT between the moment vector of the wing-tank 
combination and the angle of incidence was obtained from the relationship 

02, WT = ^ an ( m 2/ m l) WT* l ac ^ precision of determining m^ ^ and 

m 2,WT directly affects the degree of accuracy of 0 2 WT» the values of 
02 yij are no i expected to be more accurate than its components . 

The lift forces on the tank were obtained from the outputs of strain- 
gage balances, and the lift phase angles were determined by using the 
electronic counterchronograph as discussed in the section entitled 
"Instrumentation. " No correction was made to these lift measurements 
due to translation of the tank since a calculation of the force due to 
translation at the worst conditions was less than 1 percent of the mea- 
sured loads. 

The moments on the tank were obtained directly from moment balances 
and the moment phase angle was determined by using the electronic counter- 
chronograph. However, these measured moments and phase angles were the 
vector sum of the aerodynamic moments and the inertia moments. The inertia 
moments were subtracted out vectorially. In some cases, the inertia moments 
could be large compared with .the measured moments. For example, for the 
condition of |a| = 0.025 radian at M = 0, the inertia moment was 19 inch- 
pounds; whereas at M = 0.75, the inertia moment was 11 inch-pounds (at a 
lower frequency) and the measured tank moment was 25 inch-pounds. 

The lift, moment, and respective phase angles were determined for 
the wing in the presence of the tank by a vector subtraction of the tank 
forces from the forces of the wing-tank combination. 


RESULTS 


The experimental results for the measured aerodynamic forces, moments, 
and phase angles for a wing-tank combination and for the tank alone over 
the end of a low-aspect-ratio wing are given in table I. The forces, 
moments, and phase angles of the wing in the presence of the attached tank 
are included. Also, given in this table are the corresponding Mach numbers 
and reduced frequencies. In order to show trends and comparisons, the 
experimental values for the wing-tank combination and for the tank alone 
are also plotted in figures k to 11. It may be noted that an investigation 
to determine the effect of the gap as compared with no gap between the tip 
tank and the wing indicated that the gap had little effect on the measured 
forces and phase angles. 
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DISCUSSION 


The total lift and moment coefficients and their respective phase 
angles for the wing- tank configurations without a fin, with a trapezoidal- 
shaped fin, and with a delta-shaped fin are discussed first. Further 
discussion will concern the oscillating lift and moment coefficients and 
corresponding phase angles for the oscillating tank (oscillating as a 
unit with the wing in amplitude and frequency). 

It should he remembered throughout this discussion that all force 
and moment coefficients have been obtained by using a reference plan- 
form area of 1 square foot and a reference chord of one-half wing chord. 
A°constant area was chosen for the presentation of the data in order to 
facilitate the evaluation of the effects of the tank and fins on the 
oscillating forces. Thus, for example, any additional lift on the wing 
tS comSLtion due to tie additional area of the fin should appear as 
an increase in the coefficient. The actual plan-form areas were as 
follows: the exposed wing, 0.80 square foot; the tip tank, 0.69 square 

foot; either the trapezoidal- or delta-shaped fin, 0.08 square foot. 

The overall lift coefficients for the various wing-tank combinations 
are shown in figure A. Although there is some scatter in the data the 
coefficients for the delta fin seem to fall somewhat above the coeffi 
cients for the trapezoidal-fin or no-fin condition. For comparison, a 
solid-line curve representing the average values (from a .previous inves- 
tigation) of experimental lift coefficients for a wing with an aspect 
ratio of 2 (with no external store) oscillating in pitch about its mi - 
chord has been shown. Since this solid-line curve appears to be somewhat 
below the coefficients of the wing-tank combination, . it may be deduced 
that on a wing with an aspect ratio of 2, the addition of a tip tank, 
in the manner employed in this investigation, tends to increase the oscil- 
latory lift forces acting on the configuration. 

The phase angles for the wing-tank combination corresponding to the 
overall lift coefficients are shown in figure 5- It may be noted that 
the phase angles for the three conditions differ only sl ^ S J t ^; ^ e .^ Se 
angles for the condition of no fin on the tank are somewhat greater than 
for the condition of a fin on the tank. For comparison, the solid-line 
curve representing average experimental values for a wing alone with 
aspect ratio of 2 is shown in this figure. Apparently, there is a small 
decrease in phase angle on a wing with an aspect ratio of 2 when a tip 
tank is added, and the addition of a fin on the tank causes a further 
decrease in "the overall lill phase angle. 

The total wing-tank oscillating moment coefficients are shown in 
figure 6 for the three tank-fin configurations. The largest overall 
moment coefficients are obtained when the tank has no fin on it, an e 
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lowest coefficients (in these tests) are obtained with the trapezoidal- 
shaped fin attached in the range of k from 0.1 to O.5. The moment coef- 
ficients obtained when the delta fin is attached are lower than those for 
the no-fin condition. This decrease in moment coefficient due to the addi- 
tion of the fin may be caused by the moment contributed by the fin opposing 
the moment acting on the wing-tank combination. Averaged experimental 
moment-coefficient data for the wing alone with an aspect ratio of 2 are 
also shown in figure 6. A comparison of this averaged data with the 
wing— tip-tank data indicated that the addition of the wing-tip tank (with 
or without a fin) caused an increase in moment coefficient. 

The corresponding phase angles for these wing-tank moment coefficients 
are shown in figure 7. The data for all three configurations indicate some 
scatter and there is no pronounced effect of configuration on the moment 
phase angles. For comparison, a solid-line curve is shown representing 
average experimental values for a wing alone with an aspect ratio of 2. 

The moment phase angles are not materially different for the wing alone 
as compared with the wing with the tip tank. 

The forces, moments, and their respective phase angles on the tank 
while oscillating with the wing were determined for the three fin condi- 
tions (no fin, trapezoidal -shaped fin, and delta-shaped fin) and are 
forthwith discussed. The oscillating- tank lift coefficients are shown 
in figure 8. The addition of either fin increases the tank lift coef- 
ficient as might be expected when it is recalled that the coefficients 
are reduced by using the same area regardless of plan form. There is 
little variation of tank lift coefficient as a function of reduced fre- 
quency k in the range of k covered. Also shown in figure 8 as a 
solid-line curve are the results of a theory which was developed for the 
oscillating air forces on an engine nacelle. This theory (found in 
appendix IV of ref. l) accounts only for the nacelle forces acting in 
front of the quarter chord of the wing and, consequently, gives only a 
rough representation of the present configuration. A comparison of the 
experimentally determined oscillating-tank lift coefficients for the 
configuration of the tank without fins with the results of the nacelle 
theory indicates very good agreement. A comparison of the magnitude of 
the oscillating-tank lift coefficients with the steady-state conditions 
may be made by reference to a recent theoretical paper (ref. 3) which 
developed the load distributions on wings in steady flow with cylindrical 
bodies at their tips. An extrapolation of the results of reference 3 
indicates that the body at the wing tip contributes about 31 percent of 
the steady total load on the combination. From a comparison of the 
oscillating-tank lift coefficients on the tank alone (no fin) in figure 8 
with the overall coefficients in figure k, it is seen that the tank con- 
tributes approximately 27 percent of the total load in this oscillating 
case. 
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The phase angles by which the tank lift leads the tank position are 
shown in figure 9 . All phase angles increase with increasing k. The 
addition of either fin caused a decrease in phase angle. A comparison 
of the phase angles of the tank with overall wing-tank phase angles 
(fig. 5 ) indicated the angles to be about the same for the same k. 

From" a comparison of the present experimental data with the theory of 
reference 1 (nacelle theory), it is seen that the experimental phase 
angles are higher than the theoretical values for a nacelle for most of 
the k range covered. 

The moment coefficients for the tank by itself (referred to the 
midchord) are shown in figure 10. The effect of the fin shape (or pres- 
ence) is indicated, although this effect is not clearly defined because 
of the scatter in the data. A comparison of the experimental data with 
the nacelle theory of reference 1 shows that the experimental tank moment 
coefficients are considerably lower than those predicted by theory for 
a nacelle. 

The tank moment phase angles for the aforementioned moment coeffi- 
cients are shown in figure 11. The phase angles for all three conditions 
lie in a scatter band of about 12® with a mean value of -10 . The nega- 
tive phase angles between moment and position indicate a stabilizing, 
moment. From a comparison of these phase angles with those of the wing- 
tank combination (fig. 7 ) it may be seen that the tank phase angles are 
more or less constant with k, whereas the overall moment phase angles 
become increasingly negative with k. Also shown in figure 11 are. the 
results of nacelle theory. These phase angles agree fairly well with 
the measured tank moment phase angles. 


CONCLUSIONS 


The aerodynamic forces and moments acting on a wing-tank configura- 
tion, with or without fins, oscillating in pitch about the wing-root 
midchord have been measured over the Mach number range from.O.lo to O. 75 . 
Comparisons with similar data obtained in a previous investigation of the 
same wing without a tip tank indicate the following conclusions: 

1. The total lifts and moments of the wing- tank combination were 
greater than those for the wing alone. 


2. The addition of a delta- or trapezoidal-shaped fin on the tip 
tank caused a decrease in the total moments in contrast to the no-fin 
condition*, however, these moments remained greater than for the wing, 
alone. Furthermore, the moment phase angles were not changed appreciably 
by the presence of the tank, so that the addition of the tip tank with 
or without fins had a stabilizing influence on the configuration. 
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3. The overall lift and moment phase angles did not change greatly 
with the addition of the tank. 

k. The tank lift increased with the addition of either fin. 

5. The tank lift and phase angles agreed very well with the results 
of an engine -nacelle theory, whereas the moment agreement was not so good. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Langley Field, Va., August 2, 1956. 
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APPENDIX A 


CORRECTION OF ROOT REACTION FOR INERTIA AND AEROELASTIC 
EFFECTS TO DETERMINE THE TOTAL AERODYNAMIC LOAD 

In the present experiment on a cantilever wing, knowledge is desired 
of the total aerodynamic load which develops solely from the torsional 
oscillations of the wing. This aerodynamic load is not equal precisely 
to the reaction at the wing root because of the presence of secondary 
bending reactions which come as a result of the freedom of the wing to 
be excited slightly in a bending oscillation. A correction must, there- 
fore, be applied to the measured root reaction to obtain the aerodynamic 
load associated directly with the torsional motion. This appendix derives 
and shows the magnitude of this correction. The derivation is made in 
general terms of a wing of variable cross section; the correction is then 
applied to the uniform wing-tank, combination. 

On the basis of the engineering beam theory, the differential equa- 
tion for bending of the wing is 


where g is the structural damping coefficient and p is the intensity 
of the applied loading. With the choice of a strip-analysis approach, 
the loading for the case under consideration may be written 


The first term in this expression is the inertia force associated with 
the wing mass; the second and third terms refer, respectively, to the 
apparent air-mass inertia effect and the aerodynamic damping associated 
with bending oscillations; and the fourth term refers to the torsionally 
induced aerodynamic loading, which herein is regarded as the applied 
forcing function. The second and third terms were established by using 
oscillating-flow theory for two-dimensional incompressible flow as a. 
guide; the lift-curve slope a and the F and G coefficients, which 
are like the in-phase and out-of -phase Theodorsen flutter coefficients, 
are to be selected as appropriate to the case being treated. 



(Al) 


P = -my 


y - | pcV(F + iG)y + P 


(A2) 
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Since harmonic motion is involved, the loading and the deflection 
may be written, respectively. 


P - Lfxje 10 * 

(A3a) 

y = Y(x)e‘*- cu ^ 

(A3b) 


Now, with the use of equations 


(l + ig )— ir El 

dx^ dx 


(A2) and (A3), equation (Al) reduces to 
= o&aY - i pcVcu(F + iG)Y + L (a4) 


where 


m 


m + 



A convenient and fairly accurate approximate solution to this equation 
can be obtained by expressing the deflection in terms of the fundamental 
bending vibration mode of the wing; the choice here of only a single-mode 
expansion is considered adequate since the forcing frequencies used in 
the experiments were below the fundamental wing frequency. Thus, 


Y 


a l y l 


(A5) 


where a 1 represents the response amplitude to be determined and y 
is given in terms of unit tip amplitude and satisfies the equation 




2 - 


= % myi 


(a6) 


In accordance with the Galerkin procedure for solving differential equa- 
tions, equation (A5) is substituted into equation (a 4) which is then 
multiplied by y^ and integrated over the length of the wing. The result, 

with the use of equation (A6), is 

2 ITT, PT 

a- L (l + ig)%. J = a^J - ia L ^ — o^(F + iG)H + J 1^ dx (A7) 
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where 


J = / my 2 dx + idy 


h = / yi 2dx 
Jq c r 


mj. = 


rtpc T 


k = 


cue. 


2V 


m^. is the mass at the wing tip, and c r is some convenient reference 

chord which is usually taken at approximately the three-quarter-span sta- 
tion. The desired response amplitude can now be determined directly from 
equation (A7); hence, 


Iy 1 dx 


9-1 — 


K 2 - ^ 2 )j - f £ + i (g % 2 J + J l V 1 ^ 2 ) 


(A8) 


The loading on the beam is now written in terms of a^_. Substitu- 
tion of equations (A3) and (A5) into equation (A2) gives 


P = 


meje^y-^ + o> 2 y 1 - i | pcVcu(F + iG)y^ 


a^ + L/ e 


icut 


The root shear or reaction, which may be designated 7e lcut , f oun d 

by integrating this loading over the length of the beam; thus, 


ye 


icut - 


CD%i 


- i(F + 1G)| [ o 


L dx 


> e 


iajt 


(A9) 
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where 


% = / my-L dx + dl,. 

J 0 


B i = / 37 yi ^ 

J 0 r 


Substitution of equation (a 8) into equation (A9) and cancellation of the 
harmonic terms gives 


where 


y = 


x , ( C 1 + ic g> 


D 1 + 1D 2 


L dx 


J 0 


(A10) 


, h + a £ 2 c5l 

1 J Jt k J 


c 2 = - a £ 

^ jt k J 


D -3i! , aG”r H 

1 ' O) 2 * t J 


Iyq dx 


d = 


L dx 


0 


(All) 


As stated earlier, knowledge is desired of the total aerodynamic load 
that is associated with the torsional oscillations of the wing. This 
total load is found directly from equation (AlO) to be 



D n 


iD, 


C-jd. + 


7 
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This equation, when expressed in the complex notation of a modulus and 
phase angle, becomes 


where 


[ L dx 
J 0 




( C l d + D l) 


+ d 2 2 

2 + (C 2 d 4 - D 2 ) 2 



i0 5 

Ke y 


h 


tan 


-1 


1 

( C 1 d 2 - 

CgDj 

Id 

d| 

( c i D i 

+ C 2^2 j 

| +• + ^2 


(A12) 


(A13) 


Equation (A13) is the final equation sought. Thus, the magnitude of the 
torsionally induced air load is found simply by multiplying the correction 
K by the magnitude of the measured root shear; the phase angle between 
this load and the root shear is given by 0^. A word about the coeffi- 
cients C , D , and d may now be in order. All terms in these coef- 
ficients which contain the lift-curve slope a are related to the aero- 
dynamic damping effects associated with the bending oscillations. A 
comparison of the second term with the first term in D 2 , for example, 
will indicate how strong the aerodynamic damping is in relation to the 
structural damping. The nondimens ional term d may be seen to depend 
on the distribution of the air load L(x). For most practical cases, 
it is considered sufficient to evaluate this factor on the basis that 
the air load has an elliptic distribution. 


Some simplification results when the aforementioned relations are 
applied to the case treated herein, that is, to a uniform cantilever 
wing with a mass at its tip. When the assumption is made that there is 
no change in mode shape of a cantilever wing due to the mass at its tip, 
it can be shown that 


— 2 = mr , 

J = m / y n dx + irn,- = — + % 

Jo 1 4 


N]_ = m / y dx + im,. = 0.39 mT + m T 


H = f y x 2 dx = J 
J 0 4 


(Al4) 


B 


1 



dx = 0.39 t 
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Hence, 


C 1 " 


0.39 mT + 1% a G m r (0.39T) 
it k _ . Et 


mT 

— + “t 


rrup + 


c _ 

2 '« k e + % 


D , a G *r 4 

1 o£ « kai + IIlr 


4 


V 


aF 'll 


D 2 = 8 ~2~ + . H t 

co it k — + im,. 


and, for an assumed elliptic loading, 


(A15a) 


d = 0.29 (A15b) 

In order to determine the correction K and the phase angle 0^, 

equations (A12) and (A13) were used together with equations (A15) . The 
structural damping coefficient used for the wing was g = 0.008. The 
lift-curve slope was taken equal to the theoretical value of 2it multi- 
plied by the often-used aspect-ratio correction — -■ — , and was thus 

taken as it since the wing has an aspect ratio of 2. The F and G 
functions were arbitrarily chosen as those for two-dimensional incompress- 
ible flow. 

In order to gain an insight as to how aerodynamic -damping effects 
compare with structural -damping effects, the calculations were made for 
two conditions: (l) with structural damping only and (2) with both 

structural and aerodynamic damping included. No differences are noted 
in K for these two conditions. The following table shows the results 
for these two conditions for a frequency of 20 cycles per second: 


Structural damping 
only, a = 0 

Both structural and 
aerodynamic damping 

K 

03 > deg 

K 

0 5 , deg 

0.98 

0.01 

0.98 

0.05 
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Although a difference in phase angle is noted for the two damping condi- 
tions, the important item to note is that in hoth cases the phase angle 
is a negligible quantity. 
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APPENDIX B 


DETERMINATION OF THE AERODYNAMIC DAMPING-MOMENT COEFFICIENT 


In this appendix the method used in obtaining the aerodynamic 
damping-moment coefficient from the power-off decaying oscillations of 
a torsional spring- inertia system is given. By assuming a linear system, 
particularly with respect to the aerodynamic coefficients, and by using 
the concept that the structural damping moment is in phase with the 
angular velocity but proportional to the angular displacement, the dif- 
ferential equation of motion of the system may be given by 


or 


By definition, 


and 


! e a + Kg(l + ig vac )a = jtqS | + im 2 ) (Bl) 


»• C + iB 
a + a = 0 


C K s - *qs % m x 

r = r 


B _ 6y ac ~ ^ ^ m 2 


= A 


K, 


s - rt< l s g 


Equation (Bl) has a solution of the following form (see, for instance 
page 86 of ref. 4): ’ 


a = a e 


1/2 -iVi(i+A 2 ; 1/2 + "j 


(B2) 


For small values of A, equation (B2) becomes 

U e°4H 


a = a e 


(B3) 
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For the logarithmic decrement, equation (B3) a-t t - 0 becomes 


I ! O 


and, after n cycles, at t = 2itn/to, 


a 


n 


= a e 


-jtAn 


or 



it An 
e 


(B4) 


By taking the natural logarithm of both sides of equation (B4), 


A 




which is measured with the wing subjected to air flow. In the equation 
for A, 


^s^vac " * ^ 2 m 2 
K s - itqS | m 1 


(B5) 


Since 


and 


K s - itqS g m l = I e a)C 


K s _ -*-e a Vac 


then 


-m 


2 


I e of 


«qS 


et - 



(b6) 


Equation (B6) is the form used in the reduction of the damping-moment 
data in this paper. 
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TABLE I.- RESULTS OF AN EXPERIMENTAL INVESTIGATION ON A 
WING-TANK CONFIGURATION 
(a) No fin on the tank 


ro 


k 

M 

^WT 

^1,WT 

UI T 

1 — 1 

l*lw 

^1,W 

1^1 WT 

^2,WT 

I^It 

^2,T 

iMcdw 

^2,W 

Medium 

0.050 

0.69 

1.03 

-5 

0.28 

-3 

0.75 

-8 

0.72 

-3 

0.25 

-7 

0.47 

-1 


.053 

.67 

1.08 

-5 

.29 

-2 

.79 

-8 

.74 

-3 

.24 

-7 

.50 

-1 


.095 

.50 

1.02 

-2 

.26 

3 

.76 

-2 

.69 

-6 

.18 

-10 

.51 

-3 

Air 

.099 

.48 

1.12 

-2 

.28 

3 

.84 

-2 

.68 

-6 

.19 

-10 

.49 

-4 


.109 

.45 

1.05 

0 

.27 

5 

.78 

0 

.66 

-6 

.17 

-12 

.49 

-4 


.117 

.42 

1.14 

1 

.28 

6 

.86 

-1 

.66 

-7 

.18 

-13 

.48 

-5 


.127 

.60 

1.17 

-2 

.25 

2 

.92 

-3 

.64 


.27 

-7 

— 

— 

Freon- 12 

.128 

.40 

1.07 

2 

.26 

8 

.81 

1 

.63 

-8 

.14 

-10 

.49 

-7 

Air 

.143 

.58 

1.03 

0 

.23 

2 

.80 

-1 

.62 


.22 

-5 

— 

— 

Freon-12 

.151 

.55 

1.18 

3 

.25 

6 

.9 3 

2 

.66 


.28 

-7 




.159 

.33 

.98 

4 

.25 

12 

.64 

1 

.60 

-15 

.08 

-14 

.47 

-15 

Air 

.168 

.52 

I.07 

2 

.23 

8 

.84 

1 

.69 

-8 

.25 

-8 

.44 

-8 


.182 

.51 

1.06 

3 

.23 

5 

.83 

2 

.67 

-12 

.24 

-8 

.43 

-14 


.191 

.48 

I.07 

8 

.24 

9 

.83 

8 

.68 

-10 

.25 

-8 

.43 

-11 


.207 

.46 

I.07 

6 

.24 

10 

.83 

5 

• 70 

-11 

.24 

-8 

.46 

-12 


.228 

.44 

1.06 

8 

.23 

12 

.83 

7 

.67 

-11 

.24 

-8 

.43 

-13 


.254 

.41 

1.10 

8 

.24 

12 

.86 

7 

.65 

-16 

.24 

-8 

.41 

-21 


. 275 

.38 

1.06 

12 

.23 

15 

.83 

11 

.65 

-17 

.22 

-9 

.43 

-21 
*1 1 a 

Freon-12 

.316 

.34 

1.03 

13 

.24 

16 

• 79 

12 

.73 

-12 

.21 

-8 

.52 

-14 


.319 

.34 

1.03 

15 

.25 

19 

.78 

14 

.72 

-17 

.22 

-8 

.50 

-21 


.399 

.28 

1.03 

21 

.25 

24 

.78 

20 

.84 

-24 

.20 

-8 

• 65 

-29 


.458 

.26 

1.08 

23 

.25 

27 

.83 

22 

— 

-4l 

.15 

-8 

— 



.492 

.25 

1.09 

25 

.28 

29 

.81 

24 

.84 

-38 

.19 

-8 

.68 

-46 


• 555 

.21 

1.11 

30 

.30 

31 

.81 

30 

• 98 

-30 

.16 

-6 

.84 

-34 


.638 

.18 

1.20 

32 

.34 

33 

.86 

31 

1.18 

-36 

.15 

-7 

1.05 

-40 
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TABLE I.- RESULTS OF AN EXPERIMENTAL INVESTIGATION ON A 
WING-TANK CONFIGURATION - Continued 
(b) Trapezoidal-shaped fin on the tank 


k 

M 

Hwt 

0 

n.,wr 

I 1 1 T 

^1,T 

i»i» 

^1,W 

1 M ° t 'l WT 

ro 

s 

1^1 T 

^2,T 

Mv 

^2,W 

Medium 

0.066 

.074 

.084 

.096 

.122 

0.65 

.62 

.56 

.51 

.42 

1.06 

1.07 

1.14 

1.09 

1.02 

-3 

-3 

-3 

-1 

-1 

0.38 

.40 

.40 

.40 

.37 

-4 

-2 

-2 

-1 

-1 

0.68 
.67 
• 7^ 
.69 

.65 

-2 

-4 

-4 

-1 

-1 

0.60 
• 58 
•57 
.52 
.48 

-5 

-4 

-6 

-6 

-8 

0.13 

.13 

.13 

.10 

.06 

-5 

-6 

-7 

-6 

-12 

0.47 

.45 

.44 

.42 

.42 

-5 

-3 

-6 

-6 

-7 

Air 

.129 

.156 

.167 

.185 

.201 

.223 

.242 

.263 

.293 

.310 

.316 

•371 

.4l4 

.499 

.62 

• 59 

• 55 
.52 

.49 

.47 

.43 

.41 

•37 

.34 

.35 

.30 

• 27 
.23 

1.19 

1.19 

1.20 
1.16 
1.18 

1.12 

1.13 

1.13 

1.11 

1.09 

1.03 
1.08 

1.04 
1.04 

-3 

0 

-1 

1 

5 

1 

5 

5 

9 

10 

12 

14 

17 

24 

.40 

.38 

•37 

.36 

.38 

.35 

.33 

.35 

.38 

.36 

.29 

.37 

.32 

.35 

-4 

-5 

-3 

-1 

0 

0 

2 

3 

7 

10 

12 

14 

17 

23 

.79 

.81 

.83 

.80 

.80 

• 77 
.80 

• 78 

• 73 
.73 

• 74 

• 71 

• 72 
.69 

-2 

2 

-1 

1 

7 

1 

6 

6 

10 

10 

12 

14 

17 

25 

.60 

• 50 
.54 
.55 
.47 
.54 
•55 
.52 

• 53 
.55 
.62 

• 57 
.49 
.56 

-17 

-16 

-12 

-21 

-13 

-17 

-17 

-21 

-25 

-17 

-27 

-38 

-55 

.15 

.15 

.12 

.20 

•27 

.19 

.18 

.25 

.18 

.28 

.31 

.30 

.25 

.21 

-9 

-9 

-7 

-11 

-15 

-11 

-10 

-14 

-10 

-16 

-17 

-17 

-14 

-11 

.35 

.42 

.35 

.20 

.35 

.37 

• 27 

• 35 
.28 
.31 
.28 
.28 
.43 

-20 

-19 

-13 

-29 

-14 

-20 

-20 

-27 

-34 

-17 

-38 

-59 

-75 

Freon-12 


ro 
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TABLE I.- RESULTS OF AN EXPERIMENTAL INVESTIGATION ON A 


WING-TANK CONFIGURATION - Concluded 


(c) Delta-shaped fin on the tank 



I 


< t 



L-860V7.2 

Figure 1.- Photograph of wing-tank combination. 
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Figure 2.- Diagrammatic view showing various fins of wing-tank combination. 



Figure 3*~ Diagrammatic view of oscillator mechanism and wing-tank 
combination mounted in the Langley 2- by 4-foot flutter research 

tunnel. ro 
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Figure 4.- Variation of total experimental lift coefficients with 
reduced frequency for the oscillating wing-tank combination. 
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Figure 5*- Variation of total experimental lift phase angles with 
reduced frequency for the oscillating wing-tank combination. 
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Figure 6.- Variation of total experimental moment coefficients with 
reduced frequency for the oscillating wing-tank combination. 
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Figure 7 • - Variation of total experimental moment phase angles with 
reduced frequency for the oscillating wing-tank combination. 
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Figure 8.- Variation of experimental lift coefficient with reduced 
frequency for the oscillating tank. 
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Figure 9*- Variation of experimental lift phase angle with reduced 
frequency for the oscillating tank. 
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Figure 10.- Variation of experimental moment coefficient with reduced 
frequency for the oscillating tank. 
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Figure 11.- Variation of experimental moment phase angle with reduced 
frequency for the oscillating tank. 


-o 


NACA TN 5822 





